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SUMMARY 


The design and analysis of proposed supersonic cruise aircraft structures 
has required extensive use and new development of computer-aided methods. 

This paper reviews such computer-aided methods which have been and are being 
developed by Langley Research Center in-house work and by related grants and 
contracts. Synthesis methods to size structural members to meet strength and 
stiffness (flutter) requirements have been emphasized in this work and are 
described. Because of the strong interaction among the aerodynamic loads, 
structural stiffness, and member sizes of supersonic cruise aircraft struc- 
tures, these methods have been combined into systems of computer programs to 
perform design studies. The approaches used in organizing these systems to 
provide efficiency, flexibility of use in an iterative process, and ease of 
system modification are discussed. 


INTRODUCTION 


Supersonic cruise aircraft tend to be large and flexible, and realistic 
determination of their aeroelastic behavior requires finite-element structural 
analysis and sophisticated steady and unsteady aerodynamic loading analysis. 

Such analysis methods have been developed and are available in the form of 
computer programs and systems for use in supersonic cruise aircraft research; see, 
for example, references 1-4. The capability to size structural members to 
have low structural mass and retain adequate strength and stiffness to meet all 
the design requirements is also important. Computer-aided methods for struc- 
tural synthesis are receiving considerable attention, but are not yet developed 
to the same level of completeness and sophistication as computer-aided struc- 
tural analysis methods. 

The purpose of the present paper is to review computer-aided methods for 
structural analysis and synthesis which are being developed by Langley Research 
Center in-house efforts and by related grants and contracts. Methods for 
sizing structures to meet both strength and stiffness (flutter) requirements 
have been emphasized in this work and are described. Because of the strong 
interaction among the aerodynamic loads, structural stiffness, and member 
sizes of supersonic cruise aircraft structures, these methods have been combined 
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into systems of computer programs to perform design studies, Th.ese programs 
include static and dynamic structural analyses, steady and unsteady aerody- 
namics, static aeroelastic and flutter analyses, and methods for sizing 
structural members. 

The finite-element structural analysis program is the central focus in 
developing such a system. The "SPAR" computer program is currently being used 
for this purpose at Langley, and its capabilities and features are presented. 
The theoretical basis for SPAR is given in reference 5. 

Both optimality criteria and mathematical programing procedures have been 
used to size structures to meet strength requirements. These procedures 
involve Iterative processes, and the analysis and design methods have been 
tailored for efficient use in such processes. A methodology is described which 
Includes these procedures and which results in simultaneous calculation of 
aircraft design loads and structural member sizes. Efficient analysis and 
reanalysis methods and approximation techniques are used with mathematical 
programing procedures to provide capability to size large-scale structural 
models to meet flutter requirements. The initial implementation of these 
methods in a pilot computer program and subsequent incorporation into SPAR 
is discussed. 

The method of organization of these various computer programs into systems 
necessary for structural studies in supersonic cruise aircraft research is an 
important consideration in terms of time and cost required for development and 
ease and flexibility of use of the software product. Several approaches were 
used in organizing the systems described herein and their relative merits 
are discussed. 


STRUCTURAL ANALYSIS 


The emphasis in this paper is on structural synthesis or sizing. However, 
a fundamental part of a structural synthesis system is the program used for 
structural analysis. The SPAR program is currently being used for this pur- 
pose- at Langley, and its technical capabilities and features which facilitate 
its incorporation into a synthesis system are presented in this section. The 
actual use of SPAR in structural synthesis systems is discussed in subsequent 
sections . 


SPAR Technical Capabilities 

SPAR is a system of computer programs capable of computing static deflec- 
tions and stresses, natural vibration frequencies and modes, and buckling 
loads and mode shapes of linear finite-element simulations. The system is 
composed of a group of individual software "processors" which are used in a 
logical sequence to perform a desired analysis. Each processor is designed 
to perform a limited, yet distinct and complete, function. Therefore, a very 
high degree of modularity and user flexibility is provided by this sytem. 

These processors are able to communicate directly and automatically with a 
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data base complex which contains the information generated by or used by other 
processors during a computer run. A schematic of the organization of the 
SPAR system is shown in figure 1. 

The technical processors are operational on UNIVAC-1100 and CDC 6000/CYBER 
series computers. On UNIVAC systems, the processors are separate, absolute 
programs which are executed sequentially. On CDC systems, all the processors 
are designated as primary overlays in a single, absolute program which has a 
zero level overlay to call the processors for execution. This zero level 
overlay is configured to simulate UNIVAC operation; that is, the required input 
stream or file is identical to that used on UNIVAC systems. This input file, 
shown on the left-hand side of figure 1, is designed to provide for effective 
interactive operation. In interactive operation, the @XQT statements which 
calls a selected processor for execution and the related input data for that 
processor are typed in sequentially at a user console keyboard after successive 
prompts from the operating system. The computational sequence is continued by 
using another @XQT statement to call the next desired processor. 

The SPAR program is organized into two functional processor groups; tech- 
nical processors and utility processors. The functions of each of the SPAR 
processors are given in table 1. (Table 1 applies explicitly to a recent 
version of SPAR, denoted level 10.) The technical processors can be used in 
any logical sequence to perform a desired analysis such as static stress, 
vibration analysis, etc. A general characteristic of the technical processors 
is that they are very efficient with respect to both core storage and central 
processing time requirements. The method for handling the large, sparse 
matrices encountered in finite- element structural analysis in a manner to 
achieve this efficiency is described in reference 5. All of the technical 
processors make extensive use of auxiliary disk storage and operate with 
computer memory which is automatically allocated so that large structural 
simulations can be analyzed. The method used for eigenvalue and eigenvector 
calculations (ref. 6) makes it possible to perform a vibration analysis 
without first reducing the number of degrees of freedom being considered. This 
capability is useful in structural synthesis since a single simulation can be 
used for both static and dynamic calculations. Finite elements which are 
currently available in SPAR include axial bars, beams of general cross section, 
triangular and quadrilateral plates having an option to specify coupled or 
uncoupled membrane and bending stiffness, and quadrilateral shear panels. The 
provision for warping of the quadrilateral plates and an option to specify 
plate properties as layers in a laminate of composite materials are recent 
additions . 

Three types of utility processors are available in SPAR; Cl) the data 
complex utility, DCU, which is used to input, output, and manage data sets on 
the data complex, (2) the arithmetic utility system, AUS, which provides 
general matrix arithmetic capability and data set construction, and C3) on-line 
or off-line plotting. All of these utilities provide the user with an inter- 
face to the data complex. The contents and features of this data complex is 
discussed in the following section. 
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SPAR Data Complex 


The organization of the data complex provides ready access to the data 
which is produced by the various processors during an analysis run. This 
ready access of data is an important feature since it greatly facilitates the 
use of SPAR in a structural synthesis system. The data complex is composed of 
a group of up to 20 libraries, as shown at the bottom of figure 1. Data 
needed to communicate among the processors during an analysis run are stored 
on these libraries. User input to each processor specifies the library from 
which the processor will retrieve needed data and the library on which the 
processor will store the data it generates. For many analysis runs, only a 
single library is used. These libraries are direct access files which contain 
both the data and a directory to the data. Hence, they are separate, self- 
contained files which can be saved without alteration on disk storage at the 
end of an analysis run and used in subsequent runs. These files are recognized 
by the computer operating system as having names of SPARLA, . . . SPARLT, as 
indicated on figure 1. 

Each library contains data sets which are produced by the SPAR processors. 
Each data set has a four-word identifying name. For example, the set of static 
displacements for load case 10 and boundary condition set 1 would be named 
STAT DISP 10 1. Through such names, the SPAR processors are able to locate 
and access automatically all of the information needed to perform a particular 
analysis. A directory, or table of contents, is stored on each of the 
libraries in order to determine the size characteristics and location on disk 
of a data set with a specified name. An example of such a table of contents 
is shown in table 2. 

Each line in the table denotes a data set and gives the sequence it was 
written on the file, the location on disk, date and time of creation, error 
code, size and data type information, along with the name used to reference 
the data set. The data complex utility provides the capability to print these 
tables of contents, as well as individual data sets. A COPY command is pro- 
vided to transfer specified data sets between libraries. This capability is 
very useful in saving selected information between runs. The XCOPY and the 
XLOAD command are used to copy data sets from the direct access libraries to 
sequential files and vice versa. These sequential files can then be used to 
Interface other programs to the SPAR system. The arithmetic utility system 
provides for input of user defined data sets into the library, as well as 
general matrix arithmetic operation on the data sets. 

The flexibility of use of the technical processors and the ability to 
access and manage data readily in the SPAR system were very useful in 
developing structural synthesis or sizing systems. The next section of this 
paper describes two structural synthesis systems which use SPAR. 


STRUCTURAL SYNTHESIS 


Structural synthesis is the sizing (determination of areas' or thickness) 
of members in a structure to satisfy a prescribed objective, for example 
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minimum mass. Design constraints such as allowable material strength and 
flutter boundaries, which also must be satisfied by the structure, are imposed 
during this process. In this section, methods are described for the synthesis 
of structures to meet static strength requirements and flutter requirements. 

These methods were developed to perform studies of different structural 
concepts in a supersonic cruise aircraft research project. During these 
studies, the structural members are sized to meet strength requirements and 
subsequently analyzed for flutter and stiffened, if necessary. The finite- 
element representation used in these studies is shown in figure 2 to character- 
ize the level of complexity for which the methods were intended. A half model 
composed of one side of the simulation is used in the structural analysis with 
proper boundary conditions at the centerline to represent symmetric behavior 
of the complete vehicle. The half model contains 746 grid points with a total 
of 2138 unrestrained degrees of freedom and 2369 elements. Attention is 
focused on the design of the wing box structure which is modeled in greater 
detail than the remainder of the structure. During structural sizing, only 
thicknesses of members in the wing structural box (384 wing rib and spar 
shear webs and 334 wing cover panels) are calculated. Sizes of all other 
structural members are fixed at predetermined values to represent correct 
overall stiffnesses for the remainder of the aircraft. The mass of payload, 
fuel, and other items making up the total mass of the vehicle are distributed 
over the model by the analyst to provide both the correct inertia relief loads 
during maneuver and the correct mass distribution for calculating free-free 
vibration characteristics for use in flutter analysis procedures. Numerical 
results from design studies using this model are contained in reference 7. 
Herein, a description of the methods is presented. 


Synthesis for Static Strength Requirements 

Computation in strength sizing .- Supersonic cruise aircraft tend to have 
thin, large area wings and are, therefore, very flexible. Strength sizing of 
such structures poses complex computational problems since the element sizes 
depend on the stresses caused by the aerodynamic loads that are governed by 
the structural deformations which, in turn, depend on the element sizes. 

The problem is further complicated by the necessity to consider Inertia 
dependent loads such as those generated while rolling on a rough runway (taxi 
loads) and proper consideration of the jig shape (the prescribed aircraft shape 
at construction which will deform to the desired aerodynamic shape at the cruise 
condition) . The interaction of these physical quantities are indicated in 
figure 3. A computer-aided synthesis system to size a structure for strength 
inherently has the same interactions among the computer programs used to 
calculate these physical quantities. 

The two principal computational tasks that must be performed during the 
strength sizing process are (1) the calculation of aeroelastic loads, and 
(2) sizing of the structural members. These two tasks are Illustrated 
schematically in figure 4. The conventional procedure is first to calculate 
aeroelastic loads based on initial member sizes and then use these loads to 
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calculate new member sizes in a separate, sequential operation. Following the 
calculation of new member sizes, the whole procedure should be repeated because 
changes in member sizes alter the structural flexibility which impacts the 
aeroelastic loads. This iteration should continue until the changes in overall 
structural flexibility become sufficiently small to have a negligible effect on 
the aerodynamic loads. In conventional preliminary design practice, this 
iterative process is seldom fully converged because of the time and cost 
required. The interaction effect can be small if good initial member sizes 
from a previous, similar design are used in the load calculations. However, 
to reduce flow time for input preparation during preliminary design, it may be 
necessary to use crude estimates of the initial member sizes. In this case, 
large changes in structural flexibility could occur during structural sizing, 
and significant changes in the design loads could result. 

Sizing of Aerospace VEhicle Structures (SAVES) .- A system called SAVES has 
evolved from a continuing effort at the NASA Langley Research Center to develop 
structural synthesis methods to satisfy strength requirements. An initial 
version of this sytem used the conventional, sequential process of calculating 
design loads and structural member sizes (ref. 8). The SAVES system has 
been restructured into a parallel organization as shown in figure 5 so that the 
aerodynamic loads and structural member sizes are both updated at the same time. 
Simultaneous convergence of compatible design loads and structural member sizes 
is shown in reference 9 to be achieved with a substantial reduction in computer 
CPU time from the conventional procedure. 

Initially, in this process, the jig shape and maneuver deflected shape 
are assumed to be the same as the rigid cruise shape, and initial member sizes 
are specified to be constant over large areas of the structure. The process 
then proceeds in two simultaneous iterative loops as shown in figure 5. 

Loop (I), shown by double-line arrows, performs the maneuver aeroelastic load 
calculation, and loop (II), double dashed-line arrows, iteratively sizes the 
structural elements. Loop (II) uses the element sizes to generate a structural 
stiffness matrix in SPAR. Cruise, maneuver, and taxi loads are applied to 
this stiffness matrix and deflections and stresses are generated. The stresses 
are used in loop (II) by a resizing algorithm to calculate new element sizes. 

In loop (I) , the deflections from the rigid cruise shape due to cruise loads 
are used to calculate the jig shape, and the deflections due to maneuver loads 
ar^ then added to the jig shape to get the maneuver shape. The aerodynamic 
Influence coefficient matrix for the maneuver condition is saved so that 
updated loads are calculated through multiplication of this matrix by the 
panel slopes obtained from the maneuver deflections. 

There are ten computer programs currently in the SAVES system. The first 
is a preprocessor for generating finite-element structural models of the type 
shown in figure 2. A linear aerodynamic lifting surface program (ref. 2), 
the SPAR structural analysis program, and a structural sizing program are the 
primary programs used in the simultaneous calculation of design loads and 
member sizes. The sizing program uses either a weight-strength method or a 
mathematical nonlinear programing method to size individual structural 
elements, such as wing cover sandwich panels, while a fully stressed design 
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is being carried out for the overall structure. Six other programs perform 
auxiliary functions required in the process such as calculating aerodynamic 
panel slopes, updating the aerodynamic pressure distributions, trimming the 
aircraft in flight, computing the jig shape, processing the SPAR stress output, 
and preparing input data for the resizing program. Graphics programs are used 
to generate displays of design data. The SPAR structural analysis program is 
also used to calculate vibration mode shapes and frequencies for use in sub- 
sequent flutter analyses. 

The SAVES system is tailored to perform the entire strength sizing process 
with computational efficiency in a single computer run. Using this system, it 
is feasible to obtain structural design results with converged, compatible 
loads and member sizes based on sound analytical methods within the cost and 
time constraints of preliminary design. 


Synthesis for Flutter Requirements 

Structural synthesis methods to meet flutter requirements are important 
in supersonic cruise aircraft research as indicated in references 10 and 11. 

The usual design process for aircraft wing structures consists of performing a 
structural synthesis to meet strength requirements and, if required, determin- 
ing the additional stiffness or mass distribution to prevent flutter using the 
analyst's judgment in a trial-and-error procedure. In the past few years, 
there has been considerable work done to develop automated structural synthesis 
methods to satisfy flutter requirements. A comprehensive review of this work 
is contained in reference 12. In the present paper, work that. was performed in 
this area under a NASA grant (see refs. 13 to 17) will be reviewed. Initially, 
a pilot program (WIDOWAC) was used to develop, improve, and test methods that 
were needed. Currently, these methods are being combined with the SPAR pro- 
gram to provide the capacity to handle more complex structural models. 

Wi ng Design Op timization With Aeroelastic Constraints (WIDOWAC) .- The 
WIDOWAC program can be used for the synthesis of minimum mass wing structures 
which are modeled using finite-elements to represent the rib and spar webs and 
cover panels as shown in figure 6. A technique for using approximate second 
derivatives in conjunction with the Newton's method was developed for the 
optimization procedure. This mathematical programing procedure has the 
generality to allow any number of constraints to be imposed during the design 
process. Currently, multiple flutter (subsonic and supersonic) stress, strain, 
deflection, buckling, and minimum gage constraints can be included. Flutter 
modes, which exhibit a "hump" characteristic and often cause a discontinuity in 
flutter speed as the design variables are changed, are considered automatically 
in the program. 

Computation times are reduced during this process by several techniques. 
The number of degrees of the structural model are reduced by specifying that 
the wing have a symmetric cross section and modeling only either the upper or 
lower half of the structure. The finite-element structural analysis program. 
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which was developed for efficient use in this synthesis system, contains 
procedures needed for iterative analysis methods to reduce significantly 
reanalysis times compared with the original analysis of the structure. 

The number of design variables is kept small by linking the sizes of 
several structural elements under the control of a single design variable. 

Thus, the design variables do not describe the thickness of Individual finite- 
elements, but of segments of the wing, each segment containing several finite- 
elements as illustrated in figure 6. 

The wing planform is divided into a number of triangular or quadrilateral 
segments. The design variables are the thickness of the cover panels at the 
vertices of the segments. The cover panel thickness is assumed to vary 
linearly in the triangular segments. A quadrilateral segment is divided into 
two triangles, each having a linear thickness variation. Thicknesses of 
elements representing rib and spar shear webs can also be lumped under single 
design variables. 

Composite materials can be represented in the cover panels by stacking 
together orthotropic plates with stiffness properties equivalent to a fixed 
fiber orientation to construct the desired laminates. In this case, the 
design variables are the thicknesses of each of the laminae making up the 
laminate. 

WIDOWAC has been used as a test bed for evaluating optimization methods, 
constraint formulations, and analysis algorithms suitable for structural 
synthesis for. flutter. Improvements in efficiency of these synthesis pro- 
cedures and continuing decrease in computer costs have created an atmosphere 
where these methods can potentially be used on complex structural models which 
are encountered in supersonic cruise aircraft research. In the next subsection, 
a description is given of work being done to study complex models by coupling 
the large-scale, general-purpose, finite-element and data management capa- 
bilities of SPAR with the structural synthesis capabilities of WIDOWAC. 

Program for Analysis and Resizing of Structures (PARS) .- PARS is being 
developed by converting the flutter analysis and optimization capabilities in 
WIDOWAC into new SPAR technical processors. This approach will result in a 
system having the same advantages of SPAR such as modularity, flexibility of 
use, etc. , which were discussed previously. To date, the first six processors 
have been developed and are being tested (see fig. 7). The functions of these 
processors are: 

1. AERO - Defines the required aerodynamic input and the interface 
between the structural and aerodynamic grids. 

2. SUBKRN - Calculates aerodynamic matrices which are based on subsonic 
kernel function aerodynamics and are Independent of vibration modes. 

3. GAF - Calculates generalized aerodynamic forces corresponding to the 
vibration modes. 
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A, DKM - Defines the design variables for use in optimization and 
calculates the derivatives of the generalized stiffness and mass matrices 
with respect to these design variables. 

5. OPT - Performs a flutter analysis of a specified structure or per- 
forms a structural optimization to achieve a minimum mass design that satisfies 
the flutter constraints. 

6. DGRA - Calculates the derivatives of displacement vectors with 
respect to the design variables. 

The first five processors provide a flutter analysis capability for sub- 
sonic Mach numbers. A V-g flutter analysis can be performed at a given alti- 
tude or the altitude can be determined at which the flutter velocity and 
specified Mach number are compatible, or matched, in a given atmosphere. The 
flutter synthesis option currently in PARS uses the same set of vibration modes 
all through the resizing process. This requirement limits the application to 
moderate differences between the initial design and final optimized design. 

The limitation will be removed by adding capability to make it convenient to 
return to SPAR to calculate updated vibration modes after significant changes 
in the structure have occurred and then repeat the optimization in PARS. The 
last processor, DGRA, calculates the sensitivity of static displacements and 
stresses to changes in design variables. These quantities are needed to add 
displacement and stress constraints to the optimization procedure. 

The GDC version of PARS is organized in a separate program using the 
same input/output and data management routines as SPAR and communicating to a 
common data library as shown in figure 7. This arrangement is desirable during 
program development because as programing changes are made in the PARS pro- 
cessors, the entire SPAR system need not be reloaded in the computer. The 
modularity of the system has expedited the development and testing of PARS and 
forms the basis for continued developments in structural synthesis. 


ORGANIZATION OF STRUCTURAL SYNTHESIS SYSTEMS 


Structural synthesis systems are characteristically quite large because 
of the number of programs and procedures necessary to represent all the 
disciplines which interact in the sizing process. Thus, the organization of 
these large systems is an important consideration. Two types of organization 
represented by SAVES and PARS have been used for the synthesis systems pre- 
sented herein, with the principal difference being the method used to transfer 
data among programs and processors. Both of these systems are implemented on 
the GDC Network Operating System, NOS, and PARS is also implemented on 
UNIVAC EXEC 8. SAVES makes extensive use of operating system capabilities 
such as the permanent file system, procedure files, and utilities to change 
data files which are in the form of card images. 

The organization of the SAVES system is indicated on figure 8 which 
depicts a portion of the total system. The bottom level of the system is a 


collection of permanent and temporary sequential data files. Primary or basic 
input data, including finite— element structural model, aerodynamic paneling 
data, and load case information, for the system are stored on permanent disk 
files in the form of card images which are copied to local files during pro- 
gram execution. Data which are generated by one program and used by another 
such as deflections, stresses, and aerodynamic pressures are handled as tem- 
porary local files during a run. These temporary files can be saved and used 
for restart purposes in a later run. 

The next level of the system consists of decks containing file manipula- 
tion instructions which alter and combine files to form input files for the 
analysis programs. Some of the analytical programs used for calculating aero- 
elastic loads in SAVES are shown on the next level. In this organization, the 
SPAR program must interact with other separate programs and data which are not 
part of the data library. An example of the data flow which occurs in the 
series of programs shown in figure 8 will be given to illustrate the method 
of interfacing the data. 

The process is initiated using the cruise shape geometry in the program 
AEROIN to calculate angles of attack of the aerodynamic panels for use in the 
subsequent program PRESS which calculates the corresponding pressure distri- 
bution. The aerodynamic Influence coefficients which are generated in PRESS 
are saved for subsequent calculations. Equilibrated loads for the aircraft in 
flight are calculated in TRIM based on these pressure distributions. These 
loads are output in the proper form of input card images for the SPAR program, 
and "modify instruction" decks are used to merge these data into the finite- 
element model data. SPAR is then executed to calculate deflections which are 
stored on the SPAR data library. The SPAR utility, XCOPY, is then used to 
transfer these deflections to a named sequential file which can be then used 
by JIG to calculate the j Ig shape and by AEROIN to calculate a new maneuver 
shape. The process is then repeated with the data files containing the 
deflections, jig shape, and maneuver loads being updated during each iteration. 

These programs are executed in a prescribed manner by procedure files, 
which are sequences of executable control cards, as indicated by LLOOP and 
XQINIT in figure 8. A job control deck, which is the highest level procedure 
file, is used to control the overall process. Other procedure files indicated 
by the dotted lines and boxes control other portions of SAVES such as struc- 
tural sizing, vibration analysis, and subsequent flutter analysis. 

A system with the same organization as SAVES can be developed quickly and 
inexpensively on the basis of standard operating system features as indicated 
in reference 18. However, as the capabilities of a structural synthesis 
system are Increased, the number of program and data files could become quite 
large and, hence, difficult to manage and document. 

The organization of the PARS system, on the other hand, uses directly the 
data file management capabilities of the SPAR program as was indicated on 
figure 7. This bookkeeping system provides a standard form for storing or 
retrieving data sets by name using direct access files. A limited number of 
subroutines (user is primarily concerned with two) are used to read or write a 
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data set identified by a four-word name from or into a data library. The size 
(number of words), along with other descriptive information for the data sets, 
are contained in tables of contents embedded in the data libraries. Documen- 
tation Involves only specifying the physical meaning of the contents of each 
data set and specifying which of these data sets are input to or output from 
each of the processors. 

The capability for performing loops within a single input deck is cur-, 
rently not available in SPAR or PARS. Procedure files such as used in SAVES 
can be used to accomplish this needed function with the common data libraries 
being retained between multiple executions of SPAR or PARS. This organization 
provides a flexible framework on which to add improved or extended structural 
synthesis methods. 


FUTURE DEVELOPMENTS 


On-going and planned Improvements to the structural synthesis systems 
include continued work on incorporating methods to size structural components 
made of composite materials. This work will augment development of methods 
for synthesis of structures to meet overall stiffness constraints, as well as 
strength constraints. Testing of the flutter synthesis methods on large, 
complex structural models will be continued along with development of any 
needed Improvements in these methods that- may become apparent in a large 
problem environment. 

Improved calculation of the required aerodynamic quantities for these 
synthesis systems will be obtained by incorporating the capabilities of the 
SOUSSA (Steady, Oscillatory, Unsteady Subsonic and Supersonic Aerodynamics) 
program of reference 19. This program will augment or replace the steady lifting 
surface aerodynamics and unsteady subsonic kernel function aerodynamics cur- 
rently being used. SOUSSA will provide the capability to handle complex three- 
dimensional configurations with a single, common aerodynamic paneling arrange- 
ment for both subsonic and supersonic flight regimes. Some alterations of the 
program will probably be required for efficient operation in an iterative 
structural synthesis environment. This improved aerodynamic capability will 
then be organized and implemented in the form of processors for use in the 
PARS system to take advantage of the SPAR data handling and utility functions. 

The capability to assess the effects of including active control systems 
is becoming an important consideration in the development of structural anal- 
ysis and synthesis systems. This capability will be initially introduced into 
the PARS system by including the effect of active controls in the flutter 
analysis. Additional processors associated with formulating the flutter 
analysis in the proper form will be required for this work. 

As new procedures and new processors are developed for SAVES, SPAR, and 
PARS, they will provide a basic framework and flexibility of use to test and 
develop other structural synthesis methods and strategies. Therefore, the 
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consideration of modularity, efficiency, and ease of use and modifications must 
be planned and incorporated into the system so new development can concentrate 
on what is new technically and not have to rebuild or restructure existing 
programs or processors. 


CONCLUDING REMARKS 


Computer-aided design methods have been developed by incorporating sound 
theoretical procedures from a variety of disciplines into systems which have 
the capacity to handle the detailed simulations required for the analysis and 
synthesis of supersonic cruise aircraft structures. Finite-element structural 
analysis methods are available to analyze the linear behavior of a given 
structure efficiently and in great detail. Structural synthesis methods for 
sizing structures to meet strength requirements have been developed based on a 
fully stressed criterion. Compatible member sizes and aircraft design loads 
can be calculated using a simultaneous iteration procedure during this structural 
sizing process. Synthesis methods to size the members of large-scale structural 
simulations to meet flutter requirements have been developed by incorporating 
efficient analysis and reanalysis methods, approximation techniques, and design 
variable linking with mathematical programing procedures. The development of 
these systems has improved the design process by allowing more interacting 
disciplines to be considered in greater depth and with a decrease in overall 
flow time. 

Flexibility, efficiency, and ease of modification of these systems are 
important for use in advanced design applications such as supersonic cruise 
aircraft research because new demands are continually being made on these 
systems. These requirements have brought into focus the desirability of 
modular systems which make effective use of capabilities available on time- 
share operating systems. Needed technical capabilities are being developed for 
incorporation into these design systems. These capabilities include improved 
structural resizing methods (especially for stiffness considerations) and 
improved aerodynamics, as well as capability to consider composite materials 
and active controls in the structural design. 
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TABLE 1.- SPAR PROCESSORS 


Name 

TAB 

ELD 

E 

EKS 

TOPO 

K 

M 

KG 

INV 

EQNF 

SSOL 

GSF 

PSF 

EIG 

DR 

SYN 

STRP 


Function 

Creates data sets containing tables of joint locations, 
section properties, material constants, etc. 

Defines the finite elements making up the model. 

Generates sets of information for each element including 
connected joint numbers, geometrical data, material and 
section property data. 

Adds the stiffness and stress matrices for each element to 
the set of Information produced by the E processor. 

Analyzes element interconnection topology and creates data 
sets used to assemble and factor the system mass and stiffness 
matrices. 

Assembles the unconstrained system stiffness matrix in a 
sparse format. 

Assembles the unconstrained system mass matrix in sparse 
format . 

Assembles the unconstrained system initial-stress (geometric) 
stiffness matrix in a sparse format. 

Factors the assembled system matrices. 

Computes equivalent joint loading associated with thermal, 
dlslocational, and pressure loading. 

Computes displacements and reactions due to loading applied at 
the joints. 

Generates element stresses and internal loads. 

Prints the information generated by the GSF processor. 

Solves linear vibration and bifurcation buckling eigenproblems , 
Performs a dynamic response analysis. 

Produces mass and stiffness matrices for systems comprised of 
interconnected substructures. 

Computes eigenvalues and eigenvectors of substructured systems. 
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TABLE 1.- SPAR PROCESSORS CConcluded) 


Name 

AUS 

DCU 

VPRT 

PLTA 

PLTB 


Function 


Performs an array of matrix arithmetic functions and is used in 
construction, editing, and modification of data sets. 

Performs an array of data management functions including display 
of table of contents, data transfer between libraries, changing 
data set names, printing data sets, and transferring data 
between libraries and sequential files. 

Performs editing and printing of data sets which are in the form 
of vectors on the data libraries. 

Produces data sets containing plot specifications. 

Generates the graphical displays which are specified by the 
PLTA processor. 
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TABLE 2.- TABLE OF CONTENTS FOR SPAR DATA LIBRARY 

TABLE OF CONTENTS, LIBKaHY 1 

E T DATA SET NAME 


SEQ 

RR 

DATE 

TIME 

R 

WORDS 

NJ 

NI*NJ 

Y 

N1 

N2- 

N3 

N4 

1 

7 

761022 

061237 

0 

18 

1 

16 

0 

JDFl 

BTAB 

1 

8 

2 

w8 

761022 

061237 

0 

250 

250 

250 

0 

JREF 

BTAB 

2 

6 

3 

-12 

761022 

081237 

0 

12 

1 

12 

1 

altr 

BTAB 

2 

4 

4 

13 

761022 

081237 

0 

30 

2 

30 

4 

TEXT 

BTAB 

2 

1 

5 

14 

761022 

081237 

0 

10 

1 

10 

1 

MATC 

BTAB 

2 

2 

6 

15 

761022 

081237 

0 

46 

4 

46 

1 

ALTR 

BTAB 

2 

4 

7 

16 

761022 

081237 

0 

750 

250 

750 

1 

JLOC 

BTAB 

2 

5 

6 

26 

761022 

081237 

0 

250 

250 

250 

0 

JREF 

BTAB 

2 

6 

9 

32 

761022 

081237 

0 

25 

1 

25 

1 

SA 

BTAB 

2 

13 

10 

33 

761022 

061237 

0 

250 

250 

250 

0 

CON 


1 

0 

U 

37 

761022 

081237 

0 

2250 

250 

2250 

1 

OJJT 

BTAB 

2 

19 

12 

73 

761022 

081241 

0 

3136 

196 

896 

0 

DEF 

E43 

11 

4 

13 

129 

761022 

081241 

0 

2 

1 

2 

0 

GO 

E43 

11 

4 

14 

130 

761022 

081241 

0 

15 

1 

15 

0 

6TXT 

E43 

11 

4 

15 

131 

761022 

061241 

0 

12 

12 

12 

0 

nEl2 

BTAB 

1 

11 

16 

132 

761022 

061241 

0 

5 

1 

5 

0 

KE 


0 

0 

17 

133 

761022 

081241 

0 

7 

1 

7 

0 

NS 


0 

0 

18 

134 

761022 

081241 

0 

1 

1 

1 

3 

ELT8 

NAME 

0 

0 

19 

135 

761022 

061241 

0 

1 

1 

1 

0 

ELT6 

LTYP 

0 

0 

20 

136 

761022 

081241 

0 

1 

1 

1 

0 

elts 

NNOO 

0 

0 

21 

137 

761022 

081241 

0 

1 

1 

1 

0 

ELTS 

ISCT 

0 

0 

22 

138 

761022 

081241 

0 

1 

1 

1 

0 

ELTS 

NEL8 

0 

0 

23 

•139 

761022 

061241 

0 

1 

1 

1 

0 

ELTS 

LE3 

0 

0 

24 

140 

761022 

081244 

0 

1500 

250. 

1500 

• 1 

APPL 

FORC 

1 

1 

25 

164 

761022 

081244 

0 

1500 

250 

1500 

-1 

APPL 

PORC 

2 

1 

26 

166 

761022 

081244 

0 

1500 

250 

1500 

-1 

UNIT 

VEC 

1 

1 

27 

212 

761022 

081244 

0 

1500 

250 

1500 

• 1 

APPL 

MOTI 

5 

1 

26 

236 

761022 

081244 

0 

1 

1 

1 

•1 

TOT 

PORC 

5 

1 

29 

237 

761022 

081244 

0 

1 

1 

1 

-1 

CONY 

AUS 

1 

1 

30 

238 

761022 

061244 

0 

1500 

250 

1500 

-1 

OISP 

INC 

1 

1 

31 

262 

761022 

061245 

0 

6272 

250 

696 

0 

KMAP 


1067 

17 

32 

360 

761022 

061245 

0 

6960 

250 

1792 

0 

AMAP 


1675 

26 

33 

506 

761022 

081254 

0 

62720 

196 

320 

4 

E43 

EPIL 

11 

4 

34 

1466 

761022 

081248 

0 

20 

20 

20 

0 

OZR 

E43 

11 

4 

35 

1487 

761022 

081246 

0 

1500 

250 

1500 

m 1 

DEM 

DIAG 

0 

0 

36 

1511 

761022 

081302 

0 

42560 

250 

2240 

1 

K 

spar 

36 

0 

37 

2176 

761022 

061323 

0 

50176 

250 

3584 

1 

INV 

k 

1 

0 

38 

2960 

761022 

061333 

0 

1500 

250 

1500 

• 1 

STAT 

DI8P 

100 

1 

39 

2964 

761022 

081333 

0 

1500 

250 

1500 

m J 

STAT 

REAC 

2 

1 

40 

3008 

761022 

081345 

0 

10780 

196 

5555 

m 1 

STRS 

E43 

100 

1 
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INPUT FILE 
§ XQT TAB 

(TAB DATA) 

@ XQT ELD 

(ELD DATA) 

0 XQT* K 

0 XQT INV 


DATA BASE 
COMPLEX 


TECHNICAL PROCESSORS 


TAB 

TOPO 

EQNF 

EIG 

ELD 

K 

SSOL 

DR 
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M 
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PSF 

STRP 


INV 




LIBRARY 1 
DATA SET 1 
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FILE = SPARLA 


LIBRARY 2 
DATA SET 1 

= SPARLB 


UTILITY 

PROCESSORS 

AUS 


LIBRARY 20 


= SPARLT 


Figure 1.- Organization of SPAR computer program. 
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Figure 3.- Interaction of data which is calculated during 
structural sizing of flexible aircraft . 
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Figure 4.- Iterative procedures to sequentially calculate 
aeroelastlc loads and to size structural members. 





















Figure 5.- Schematic of procedure for simultaneous 
calculation of aircraft design loads and struc- 
tural member sizes. 
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Figure 6.- Characteristics of WIDOWAC flutter synthesis 

computer program. 








Figure 7.- Schematic showing processors contained in 
PARS computer program and its interaction with 
SPAR program. 



Figure 8.- Organization of a portion of 
SAVES structural sizing system. 


















